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The dark matter enclosed in a density perturbation with a large initial amplitude (5p/p > 10 -3 ) 
collapses shortly after recombination and forms an ultracompact minihalo (UCMH). Their high 
central densities make UCMHs especially suitable for detection via astrometric microlensing: as 
the UCMH moves, it changes the apparent position of background stars. A UCMH with a mass 
larger than a few solar masses can produce a distinctive astrometric microlensing signal that is 
detectable by the space astrometry mission Gaia. If Gaia does not detect gravitational lensing by 
any UCMHs, then it establishes an upper limit on their abundance and constrains the amplitude of 
the primordial power spectrum for fc — 3500 Mpc -1 . These constraints complement the upper bound 
on the amplitude of the primordial power spectrum derived from limits on gamma-ray emission 
from UCMHs because the astrometric microlensing signal produced by an UCMH is maximized if 
the dark-matter annihilation rate is too low to affect the UCMH's density profile. If dark matter 
annihilation within UCMHs is not detectable, a search for UCMHs by Gaia could constrain the 
amplitude of the primordial power spectrum to be less than 10 -5 ; this bound is three orders of 
magnitude stronger than the bound derived from the absence of primordial black holes. 



I. INTRODUCTION 

Structure formation is hierarchical: small-scale density 
perturbations form the first dark matter minihalos, and 
these minihalos are later absorbed into larger dark mat- 
ter structures. If these small-scale density perturbations 
have the same average initial amplitude (Sp/p ~ 1CU 5 ) as 
the large-scale density perturbations probed by the cos- 
mic microwave background (CMB) and observations of 
large-scale structure, then the first dark matter minihalos 
form long after the Universe became matter-dominated, 
at redshifts z ^ 60 0, 51- Density perturbations with 
larger initial amplitudes can form dark matter halos 
much earlier; if (5p/p ~ 10~ 3 ) when the perturbation en- 
ters the Hubble horizon, then the dark matter within the 
overdense region will collapse to form a minihalo shortly 
after recombination (z ~ 1000) [3(, and larger fluctua- 
tions can form minihalos while the Universe is still ra- 
diation dominated |^|. Since these minihalos form in a 
denser environment than later-forming minihalos, they 
have high central densities, and they have been dubbed 
ultracompact minihalos (UCMHs). 

The abundance of UCMHs measures the primordial 
power spectrum of density fluctuations because UCMHs 
form from density perturbations with large initial ampli- 
tudes. The upper limit on the amplitude of the primor- 
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dial power spectrum derived from the absence of UCMHs 
is extremely valuable because UCMHs probe the primor- 
dial power spectrum on scales far smaller than those 
accessible via observations of the CMB and large-scale 
structure. As detailed below, UCMHs offer substantial 
improvements over other bounds on the amplitude of the 
small-scale power spectrum, but all existing limits on 
the abundance of UCMHs are derived from their poten- 
tial emission of dark- matter annihilation products 
In this paper we present a new method of searching for 
UCMHs that can provide strong constraints on the am- 
plitude of the primordial power spectrum even if dark 
matter does not self-annihilate. 

On cosmological scales, the CMB and large-scale struc- 
ture provide a direct probe of the primordial power spec- 
trum. The amplitude of temperature fluctuations in the 
CMB is proportional to the amplitude of the primordial 
curvature power spectrum on wavelengths of ~10 Mpc 
to -20 Mpc [I-Qjl. On these scales, the CMB indicates 
that the primordial power spectrum of curvature per- 
turbations is nearly scale-invariant with an amplitude 
of — 2 x 10~ 9 [3, [HI- Measurements of the matter 
power spectrum inferred from observations of lar ge-s cale 
structure (nl - fl9| and weak gravitational lensing |2"oT - [24j 
probe similar scales as the CMB, while observations of 
the Lyman-cv forest reach slightly smaller wavelengths 
(—2 Mpc) [25]. All these observations are also consis- 
tent with a nearly scale-invariant spectrum of primordial 
fluctuations with an amplitude of— 2xl0~ 9 [IE EE 53- 

These measurements of the primordial pow er spectrum 
are often cited as evidence for inflation |27H29| because 
many inflationary scenarios predict that the primordial 
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power spectrum should be nearly scale- invariant over a 
wide range of scales (see Ref. [30| for a review). How- 
ever, the power spectrum deviates from scale-invariance 
in several inflationary models. Features in the inflation- 
ary potential, such as a steps, kinks, or bumps, enhance 
the power spectrum at specific scales [3ll - l35l |. Multifield 
inflationary models can also produce steps and other fea- 
tures in the power spectrum 36-39]. Interactions be- 
tween the inflaton and other fields generate deviations 
from scale-invariance (40j ; particle production durin g in - 
flation produces a bump in the power spectrum [4ll443j . 
and excess power on small scales is generated if the infla- 
tion has an axial coupling to a g aug e field 0, 5H ■ Sev- 
eral other inflationary models [46j-|50| , including running- 
mass inflation [5lT[53l | and hybrid inflation (54T[56l |. also 
predict enhanced perturbations on small scales. Small- 
scale perturbations can also be significantly amplified af- 
ter inflation by the QCD phase transition [57j or a non- 
standard thermal history [58| . 

Clearly, it is essential that we measure the amplitude of 
the power spectrum on small scales. Other than UCMHs, 
the only probe of density perturbations on wavelengths 
smaller than ~1 Mpc are primordial black holes (PBHs). 
PBHs form when density fluctuations with an initial am- 
plitude of Sp/p <; 0.3 enter the horizon [59j, |60j. The 
abundance of PBHs is tightly constrained over a wide 
range of PBH masses; these constraints imply an up- 
per bound on the amplitude of the primordial curvature 
power spectrum of 0.01-0.06 for wavelengths between 600 
Mpc and 10~ 16 pc [61]. Limits on the abundance of 
UCMHs can provide far more powerful constraints on 
the amplitude of the primordial power spectrum because 
less extreme density perturbations are required to form 
UCMHs than are necessary to form PBHs. 

If the dark matter is thermally produced in the early 
universe, then it self-annihilates [62|-[64j], and the high 
density of dark matter within UCMHs enhances the 
annihilation rate. It was quickly recognized that the 
Large Area Telescope on the Fermi Gamma-Ray Space 
Telescope (Fermi-LAT) [65| could detect UCMHs as 
gamma-ray sources Q and that measurements of the dif- 
fuse gamma -ray background constrain the abundance of 
UCMHs H m, H3- Since UCMHs form shortly after 
recombination, gamma-ray emission from UCMHs can 
also have a profound effect on the ionization history of 
the Universe @, H3, [6{| . Consequently, measurements of 
the optical depth to the surface of last scattering also 
limit the fraction of dark matter that may be contained 
in UCMHs. Finally, the fact that Fermi-LAT has not 
yet detected gamma-rays from dark matter annihilation 
within an UCMH puts a strong upper bound on their 
number density within the Milky Way, which implies an 
upper bound on the amplitude of the primordial power 
spectrum; this analysis was initially proposed by Ref. [1] 
and was recently refined and extended by Ref. @. If 
dark matter self-annihilates, then these limits lower the 
upper bound on the primordial power spectrum by sev- 
eral orders of magnitude; if the mass of the dark matter 



particle is less than 5 TeV and it is a standard thermal 
relic, then UCMH abundance constraints imply that the 
amplitude of the primordial curvature power spectrum 
does not exceed 2 x 10~ 6 - 2 x 10~ 7 on wavelengths be- 
tween 1.6 Mpc and 0.3 pc [8(. 

Unfortunately, the numerous searches for dark mat- 
ter annihilation products have not yet detected any ev- 
idence that dark matter self-annihilates [7fj| - l78j . There 
may be an asymmetry between dark matter particles and 
antiparticles (see e.g. [79M82| ) . or the cross section for 
dark matter annihilations may be much lower than the 
value expected for a thermal relic in the standard scenario 
83, 84]. The only guaranteed signatures of UCMHs are 
gravitational; to obtain mo del- independent constraints 
on the amplitude of primordial power spectrum, we must 
search for UCMHs via their gravitational effects. 

UCMHs can be detected through photometric mi- 
crolensing Q. When an UCMH passes in front of a star, 
it produces a light curve that is similar to the light curve 
produced by a compact object (point lens) with a mass 
equal to the mass enclosed in the UCMH's Einstein ra- 
dius (assuming that the core radius of the UCMH is much 
smaller than its Einstein radius) . Currently, photometric 
microlensing searches tell us that point-like objects with 
masses between 10 -6 M Q and 10 Mq contain less than 
3% - 10% of the dark matter within the Milky Way halo 
[85H87j . However, despite its compactness, only a small 
fraction of a UCMH's mass is enclosed within its Ein- 
stein radius. If the lensed source is located in the Large 
Magellanic Cloud, a UCMH with less than 10 M within 
its Einstein radius has more than 98% of its mass outside 
this radius at its formation, and this percentage increases 
as the UCMH accretes more matter. Therefore, photo- 
metric microlensing searches do not limit the abundance 
of UCMHs because these constraints are satisfied even 
if all dark matter immediately collapses into UCMHs at 
the time of matter-radiation equality. 

In this paper, we consider astrometric microlensing by 
UCMHs, and we evaluate the constraints on the primor- 
dial power spectrum that could be obtained from a search 
for UCMHs using high-precision astrometry. When an 
UCMH passes in front of a star, the location of that star's 
image will move as the angular separation between it 
and the UCMH changes. The image trajectory produced 
by a diffuse lens is easily distinguished from the astro- 
metric microlensing signature of a point mass, offering a 
distinctive way to detect dark matter substructures. In 
Ref. [HI , Erickcek and Law explored possibility of detect- 
ing astrometric microlensing by conventional dark mat- 
ter subhalos; unfortunately, the number density of dark 
matter subhalos predicted by numerical simulations of 
galaxy-sized dark matter halos is too small for a blind 
search for astrometric microlensing events by subhalos 
to be successful. However, Ref. [88[ noted that the prob- 
ability of detecting astrometric microlensing by subhalos 
is greatly enhanced if the subhalos are more compact and 
more numerous than predicted by numerical simulations. 
Since UCMHs may be abundant in our galaxy and have 
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higher central densities and steeper density profiles than 
typical, later-forming minihalos [3j, they are optimal tar- 
gets for detection by astrometric microlensing. Never- 
theless, the image deflections produced by UCMHs are 
measured in microarcseconds, and high-precision astrom- 
etry is required to detect their astrometric microlensing 
signatures. 

Recent years have seen great growth in high-precision 
astrometric measurement capabilities using both ground 
and space-based instruments. Current all-sky astro- 
metric catalogs derived from ground-based observations 
achieve tens-of-milliacrsecond precisions (e.g. [89l . l90l|). 
while from space the Hipparcos satellite [91| provided 
milliarcsecond-precision astrometry for stars brighter 
than 9th magnitude across the entire sky. 1 For small 
fields, ground-based and space-based long-term astro- 
metric monitoring campaigns can achieve milliarcsecond- 
precision astrometry for much fainter objects (e.g. (93l — 
l98j). The development of adaptive optics systems has 
boosted the optical ground-based astrometric precision 
still further for some targets; large telescopes can ac hieve 
100-microarcsecond precision on single targets [9^, Il00j | 
while new laser-guide-star systems designed f or smalle r 
telescopes (in particular, the Robo-AO system jl0l| - |l03j ) 
offer the opportunity to cover large numbers of targets 
with similar precisions. Multi-conjugate AO systems for 
large telescopes offer the pos sibil i ty o f this performance 
over much larger fields (e.g. jl04| - |l07| V Optical and ra- 
dio interferometry have demonstr ated still higher astro- 
metric precisions (e.g. |108l4ll0l |) but these techniques 
generally require optically or radio-bright sources. 

Although these capabilities would be very useful for 
follow-up, a sensitive search for UCMHs requires all- 
sky coverage at the tens-of-microarcsecond level. For- 
tunately, the Gaia missio n, pl anned for launch in 2013, 



offers these capabilities The satellite contains 



two 1. 45m x 0.5m mirrors imaged onto a common focal 
plane containing a 1-GigaPixel astrometric camera, pre- 
cision spectrophotometers and a radial velocity spectro- 
graph. The instrument is designed to achieve end-of- 
mission astrometric accuracies of 5-14 microarcseconds 
for stars brighter than 12 th magnitude (^25 microarcsec- 
onds for 15 th magnitude stars). The important metric for 
astrometric transient searches is the single-epoch mea- 
surement precision; for Gaia, this corresponds to a sky- 
averaged measureme nt pr ecision of 23-microarcseconds 
for the ~5 x 10 6 stars |ll2j brighter than 12 th magnitude 
in the Gaia passband. 2 

Inspired by the Gaia mission, we evaluate the con- 
straints on the UCMH abundance that could be obtained 
from monthly observations of 5 million stars over six 
years, and we use these constraints to forecast bounds 
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on the amplitude of the primordial power spectrum. We 
begin in Section|TT]by reviewing the properties of UCMHs 
with particular attention to their density profiles. In Sec- 
tion IIII1 we present the image trajectories produced by 
astrometric microlensing by UCMHs, and we show how 
these image trajectories depend on the UCMH density 
profile. We then describe our detection strategy and de- 
fine lensing cross sections: the area of the sky surround- 
ing an UCMH in which a star would be detectably lensed. 
We use these lensing cross sections to calculate the proba- 
bility that surveys with varying levels of astrometric pre- 
cision would detect astrometric microlensing by UCMHs. 
In Section ITVl we calculate the constraints on the UCMH 
abundance that could be obtained by high-precision as- 
trometric surveys, and we translate these constraints to 
upper bounds on the amplitude of the primordial power 
spectrum on small scales. These bounds complement the 
bounds obtained from gamma-ray searches for UCMHs, 
and this connection is explored in Section [V] Finally, we 
summarize and discuss our results in Section [VII 



II. PROPERTIES OF UCMHS 

UCMHs form when over-dense regions with 
Sp/po <; 10~ 3 enter the Hubble horizon. The ini- 
tial mass of a UCMH, Mi, equals the entire dark matter 
content of this over-dense region at horizon entry; if R 
is the comoving radius of the overdensity and a is the 
scale factor, then the over-density enters the horizon 
when aR — H , and 
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where Hq = 100ft. km s _1 Mpc -1 is the Hubble constant, 
and O x the present-day ratio of the dark matter density 
to the critical density. Throughout this paper, the sub- 
script "hor" indicates that the quantity is to be evaluated 
at horizon entry. 

UCMHs do not grow significantly prior to the redshift 
of matter-radiation equality (z cq ~ 3250) jl!3| |. After 
matter-radiation equality, an isolated UCMH grows lin- 
early with the scale factor: 



M h {z) = Mi 
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Previous work [5j, [6( has adopted z = 10 as the redshift 
after which hierarchical structure formation prevents fur- 
ther accretion, leading to a UCMH mass growth factor 
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_ M h (z = 0) 
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300. 
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However, this much accretion is only possible if there is 
enough free dark matter to be accreted for each mini- 
halo to grow by g. If a fraction / eq of the dark matter is 
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contained in UCMHs at matter-radiation equality, then 
icq & 1/300 is required for g = 300. For larger values of 
/ C q, g is limited by the availability of free dark matter; 
in general, g = min{/ eq _1 , (1 + z oq )/ll ~ 300}. We note 
that, if UCMHs are sufficiently numerous to accrete all 
the dark matter at high redshift, it is likely that UCMHs 
will interact with other UCMHs before they are absorbed 
in larger dark matter halos. We assume that these inter- 
actions do not affect the central region of the UCMHs' 
density profiles; we revisit this assumption in Section IPV1 

The theory of self-similar secondary infall [3j, Ill4 1 1 1 5| 
models the radial accretion of collisionless particles from 
a uniform background onto a point mass and predicts a 
power- law UCMH density profile p(r) oc r -9 / 4 . How- 
ever, the angular momentum of infalling dark matter 
particles in any realistic cosmological context will be non- 
negligible jl!6l | and will lead to a significantly shallower 
profile in the central regions of these minihalos. Dark 
matter annihilation can also limit the maximum central 
density of UCMHs. Since the dynamic astrometric mi- 
crolensing signal is strongly dependent on the radial in- 
dex of the density profile [88j . we conservatively assume 
that both of these mechanisms lead to a constant density 
core with radius r c . 

Wc first consider the case of a core resulting from non- 
radial infall. Assuming that f R „ < 1/300 , the truncation 
radius of a UCMH is given by [ill, El 
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In the absence of a core, the density profile of a UCMH 
with initial mass Mi is 



p(r) 
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The density profile is time-independent because the ma- 
terial accreted after the UCMH forms has a high radial 
velocity when it reaches the inner regions of the halo and 
does not contribute sig nifica ntly to the time-averaged 
density in those regions 
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When the UCMH forms, the dark matter has a small 
velocity dispersion; thus the particles do not fall into the 
UCMH on perfectly radial orbits. If we assume that the 
dark matter velocity dispersion at UCMH formation is 
negligibly affected by the presence of the UCMHs them- 
selves, then an infalling dark matter particle's initial tan- 
gential ve locity at a radius r t from the center of the 

ucmh is fna 



<j ± ~ - 1.9 x 10~ 

77 



1000 



-0.78 / j^j- \ n.2.s 
Mr 



km; 



The radial infall approximation breaks down when the 
tangential velocity of a particle falling from r t exceeds the 



local Keplerian orbital velocity. Using the density profile 
in Eqn. ([5]) and conservation of angular momentum, this 
occurs at a radius 

1/7 

(8) 
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We take r Cinr to be the radius of the constant density core 
caused by non-radial infall. The radius of the non-radial 
infall core can be expressed solely in terms of the redshift 
of evaluation z„: 
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This radius is usually evaluated at z a of 1000, which cor- 
responds to the redshift at which perturbations of am- 
plitude 10~ 3 have grown to reach the critical overdensity 
for collapse and the UCMH forms 

We now consider cores produced by the annihilation of 
dark matter. The maximum central density of an ultra- 
compact minihalo of age t is 
m x 

Pmax 
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pc. (5) where m x is the mass of the dark matter particle and 
(av) is the product of the annihilation cross section and 
relative DM particle velocity |5| . The radius of the anni- 
hilation core is equal to the radius at which the UCMH 
density reaches p max ; from Eq. ([SJ we see that 

, / Mi \ 1/3 / m v \ - 4 /9 

r c ann = 3.3 X 10~ 4 I £ X 

\M. V 1 00 GoY i 
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pc. 
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The prevailing core radius in the UCMH density profile 
is simply the greater of r Cjllr and r Ciann . The final cored 
UCMH density profile we consider for microlensing is 

/ \ -9/4 

P(r) = Po f 1 + - J • (12) 

We set the density to zero for r > r t . If dark matter 
annihilation is the prevailing factor in determining the 
core radius (r c ann > r c nr ) then p is simply the max- 
imum central density of the UCMH p max . If the core 
due to non-radial infall prevails, then we first determine 
the overall normalization factor po,m- such that the the 
mass of the halo at matter-radiation equality Mi is equal 
to the volume integral of this density profile within the 
truncation radius at z cq , assuming that there is no an- 
nihilation. In some cases, /?o,nr ma y be greater than the 
maximum density allowed by annihilation. So in the case 
where r c nr > r c ann , we take po to be the lesser of po,nr 
and p max . 
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III. ASTROMETRIC MICROLENSING BY 
UCMHS 

A. Image Trajectories 

In gravitational microlensing by a spherically symmet- 
ric thin lens, the apparent image of a background star is 
deflected from its true position by an angle 



4s 



c 2 i 



(13) 



where d\ s is the distance between the lens and source 
planes, d s is the distance between the observer and source 
planes, £ is the impact parameter of the light ray in the 
lens plane with respect to the center of the lens, and 
£ = £/£ points from the lens center to the star. M2d(£) 
is the mass enclosed by a cylinder of radius £ collinear to 
the optical axis. Using Eq. (fT2|) . it is given by 



M 2D (£) = 4tt 
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Since we are only interested in lensing events in the 
Galactic halo, we can take d\ a /d s = 1 — (di/d a ), where 
di is the distance from the observer to the lens plane. 
In the weak microlensing regime, the angular separation 
between the true image of the star and the center of the 
lens /3 is much greater than the deflection angle a. Since 
these two angles are always collinear, we can approximate 
£ = di(a + f3) ~ d[/3. In this case, Eqn. (|13l) becomes a 
simple function for a(/3). 

Figure [1] shows the trajectories of the apparent images 
of a star at 5 kpc for varying core sizes and impact pa- 
rameters. We choose the co-ordinate system such that 
the x-axis is parallel to the proper motion of the UCMH, 
and the true position of the star is at the origin. Here, the 
separation angle /3 can be decomposed into an impact pa- 
rameter (3 y that is constant throughout the lensing event, 
and a perpendicular component j3 x along which all the 
lens motion during the event occurs. Four years of the 
lensing event are shown in Fig. [1] with points are plot- 
ted every 27 days. In the orientation shown, the UCMH 
starts at the left and moves right, while the image of the 
star starts right of its true position and eventually moves 
left. When (3 is still large, the image motion of the star 
is slow; the image accelerates until it reaches its maxi- 
mum velocity when the UCMH is directly below the star, 
precisely when when the magnitude j3 is minimized with 
P x = 0. 

It is easily seen from Eqn. (TTBl that more massive 
haloes generate a larger overall lensing trajectory. How- 
ever, the presence of a larger constant density core in such 
haloes also changes the trajectories' shapes compared to 
those of their smaller, cuspier counterparts. As shown in 
the upper plot of Fig[TJ the magnitude of the astromet- 
ric deflection a, its velocity, and its acceleration are all 
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FIG. 1: Astrometric microlensing trajectories for a star at 
5 kpc and a UCMH at 50 pc with a tangential velocity of 
200 km s -1 . The initial mass of the UCMH is M 4 = 11M . 
The true position of the star is at the origin and the UCMH 
is moving horizontally below the star; as the UCMH moves 
from left to right, the apparent image of the star moves from 
right to left. The four years surrounding the moment of clos- 
est approach are shown, with points plotted every 27 days. 
The upper plot shows the effect of varying core size while 
fixing the impact parameter at 0.4 arcseconds. The highest 
trajectory corresponds to a core caused by the angular mo- 
mentum of infalling material only, while the lower three cor- 
respond to annihilation cores with m x = 100 GeV and (a v ) = 
{3 x 10" 30 ,3 x 10~ 28 ,3 x 10" 26 }cm 3 s _1 . These correspond 
to core sizes {2.8 x 10~ 6 , 1.2 x 10~ 5 ,9.3 x 10" 5 ,7.2 x 10~ 4 } 
parsecs. The lower plot shows the effect of varying impact 
parameter with a fixed angular momentum core for the same 
UCMH configuration. In descending order, the trajectories 
correspond to impact parameters of j3 y = {0.1,0.4,1,4,10} 
arcseconds. 



diminished as r c increases if the halo mass is fixed. The 
presence of a core dramatically reduces the lensing mass 
M2D when £ <5 r c , but the relative effect of the core on 
M2D is smaller for larger £. Thus the vertical extent of 
the trajectory, which is determined when £ is small, is 
more highly suppressed than its horizontal extent, which 
is determined when £ is large. This leads to a flattened 
lensing trajectory compared to that of a smaller core. 

Given a particular UCMH, large impact parameters 
P y pose two challenges for astrometric microlensing de- 
tection. First, increasing the impact parameter beyond 
r c entails an overall reduction in the amplitude of the im- 
age's vertical deflection. Second, since we do not know 
the true position of the star, any detection scheme is 
highly dependent on the velocity of the apparent image, 
which is also reduced as j3 y increases. This is shown in the 
lower plot of FigfTJ where the size of the core is fixed but 
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the impact parameter varies. As /3 y increases, the deflec- 
tion angle is reduced (for f3 y di > r c ) and image motion 
slows dramatically - similar to the effect of increasing the 
core size. Consequently, any detection scheme for lens- 
ing events is strongly biased to those with small impact 
parameters. 



B. Detection Strategy and Lensing Cross Section 

We use the same detection technique as in Ref. [88j . 
We assume a 6 year mission lifetime with 82 epochs per 
target: a typical Gaia observing scenario. As the true 
position of the lensed star is unknown, we must search 
for anomalies in its proper motion. The first two years 
are used as a calibration period, after which we fit the ob- 
served star positions for the effects of parallax and proper 
motion. The lensing signal S is then simply the differ- 
ence between the star's extrapolated position from the 
fit and the lensed position: 
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where iV cpoc h s is the number of epochs and (Xij, Y;^) 
and (-Xp,i, Y Pi i) are the lensed position of the star and the 
predicted position of the star at each epoch, respectively. 

To calculate the signal, we establish a Cartesian co- 
ordinate system on the sky with normalized impact pa- 
rameter j3 = f3ydi/ (vt iobs) on one axis and phase, defined 
ip = /3 x ,odi/(v t t ohs ), on the other. Here, /3 Xi0 is the /3 X 
co-ordinate of the UCMH at the start of the observa- 
tional period, Vt is the transverse halo velocity, and i bs 
is the length of the observational period. Throughout our 
analysis, we take v t — 200 km s _1 and t ^ s — 4 years. The 
significance of the phase is that it denotes the proportion 
of the observational period that passes before the angular 
separation between the star and lens reaches its minimum 
(f3 x = 0). Each point on the j3 — ip plane corresponds to 
a particular initial star-lens geometric configuration for 
which we calculate the lensing signal S as defined above. 

The lensing cross section for a single UCMH Aiensed 
is the area on the sky surrounding the halo within 
which a lensed star's signal would exceed the statisti- 
cally significant minimum for detection S m i n , determined 
by the astrometric instrument employed. From numer- 
ical simulations of the detection technique |88j| . S m i n = 
1.47 SNR<Ti ns t, where SNR is the desired signal-to- noise 
ratio and o~i ns t is the single-epoch instrumental astromet- 
ric uncertainty. Given a particular UCMH and S m i n , it is 
a simple matter to determine the area on the j3 — ip plane 
that generates S > S m - ln and convert this to a solid angle 
on the sky Ai onsc( j. Since a cx (1 — di/d s ), f3 x cx ip/di and 
j3 y oc fi/di, we have Abused cx (1 - d t /d s )d^ 2 . 

We also consider a reverse calibration method that 
could be implemented given a complete 6-year dataset 
from an astrometric survey. The time-reversed motion 



of the stars during last two years are used to calibrate 
a reverse proper motion fit, which is then used to detect 
lensing events during the first four years. This process al- 
lows us to detect lensing events whose signals are largely 
generated during initial calibration period. In practice, 
the contribution to the lensing cross section from lensing 
events with phases 0.5 < ip < 1 is effectively doubled, 
neatly accounting for those lensing events for which the 
lens passes the star during the first two years of observa- 
tion. 

We restrict our forward lensing cross section to only 
include stars that reach their point of closest approach 
to the UCMH's center during the 4 year observational 
period (0 < <p < 1), allowing us to distinguish UCMH 
microlensing events from sources on binary orbits by the 
distinctive acceleration patterns in their astrometric tra- 
jectories. We note that a more sophisticated detection 
technique could potentially be subject to a weaker restric- 
tion. Furthermore, if j3 y <; r t then M2d{£) is constant, 
making the lensing event undistinguishablc from that in- 
duced by a point mass. Therefore we also restrict our 
UCMH lensing cross section within a circle of radius r t 
in the lens plane. In practice, however, the cross-sections 
for the values of S m i n and M, considered always fall well 
inside of r t . 

Individual UCMH lensing cross-sections are shown in 
Fig. [2J Given (av)/m x and S m - ln , there is a minimum 
UCMH mass that generates the required signal to be de- 
tected. The area increases rapidly as Mi increases beyond 
this minimum before settling into a weaker power-law de- 
pendence. Fixing S'miji) an increase in the core size leads 
to a greater required halo mass for detection. Finally, a 
reduction in S m - ln increases the lensed area of any given 
halo. 



C. False positives 

The lensed stars' trajectories are near-ellipsoidal on 
very long timescales, with the stars returning to their 
original positions (within ~10/ias) after a few hundred 
years (for Mi = 5M ) to tens of thousands of years (for 
Mi — 1000 Mo). In our relatively short six- year obser- 
vation window we can only see a segment of the total 
trajectory. The fastest-moving part of image's trajectory 
necessarily occurs during the observation period; away 
from that time the motion is very small (see Fig. [lj . Con- 
sequently, the trajectories induced by UCMHs are very 
different from the near-circular astrometric trajectories 
induced by point-source microlensing; these trajectories 
return to the star's tr ue p osition much more quickly. (We 
refer readers to Ref. [88( for a detailed discussion of this 
and other possible false positives.) Because of its range 
of possible trajectories, orbital motion is the most im- 
portant false positive for astrometric motion produced 
by UCMHs. However, most orbital motion can be im- 
mediately distinguished by either its repetition (for the 
shorter periods), or by the production of an anomalous 
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FIG. 2: Lensing cross sections in square arcseconds as a function of initial UCMH mass Mi with di — 50 pc and d s — 2kpc. 
The four curves correspond to different values of S m in- The three panels show the constraints for three different UCMH density 
profiles. In the left panel, the core radius of the UCMH is given by the breakdown of radial infall, while in the right two panels, 
the core radius is determined by dark matter annihilation. The dark matter particle mass is assumed to be 100 GeV, and the 
core radius is determined by the ratio {av)/m x . 



acceleration during the calibration period. 

We conducted extensive Monte- Carlo simulations to 
explore the types of orbits that could mimic the lensing 
trajectories. If the detection strategy rejects sources be- 
cause of an anomalous acceleration during the calibration 
period, we found that no Keplerian orbit was capable of 
producing a false positive signal (for eccentricities <0.99; 
see Ref. [88| for details). Longer-period highly eccentric 
orbits were rejected because of a detectable signal during 
the calibration period. Shorter-period highly eccentric 
orbits could produce a small enough acceleration to es- 
cape notice during calibration, but had orbital periods 
that would lead to a repeated motion within the Gaia 
mission. 

All microlensing-produced trajectories also produce 
some motion during the calibration period, but only the 
events in which the star and the UCMH reach their clos- 
est point early in the observational period (i.e. events 
with small phases) produce calibration-period accelera- 
tions that are detectable by Gaia. The threshold for 
calibration-period rejection of a target should be set after 
an optimization of the number of rejected false-positives 
compared to ignored microlensing events. As part of an 
improved detection strategy, which we leave for future 
work, this optimization could include further methods of 
false-positive rejection such as radial velocities, proba- 
bilistic assessments based on the proper motion and par- 
allax of the target, and further ground-based astrometric 
monitoring of the most promising cases. 

Keplerian orbit false-positives could also be distin- 
guished by Gaia directly, without the astrometric cali- 
bration period. Keplerian orbital motion, likely produced 
by a central, bright star being moved by a distant fainter 
companion, has attendant radial velocity variations. Us- 
ing the distance to the target star, a typical value for its 
mass, the size of its apparent astrometric motion, and 
the apparent period of the motion for a typical UCMH- 
like signal, we can estimate the size and orbital radius 



of a companion star that is necessary to produce the as- 
trometric motion. In almost all cases we find that the 
companion would have to be a low-mass brown dwarf. 
Using the estimated masses and orbital periods, we find 
that all astrometric orbits that could mimic a lensing 
trajectory produce radial velocities changes in excess of 
0.75 km s _1 on Gaia mission lifetime timescales. Even 
without an astrometric calibration period, many of pos- 
sible false-positive orbits can thus be directly rejected by 
Gaia's Radial Velocity Spectrometer, which will achieve 
km/s precision on these bright targets. 

Finally, a detected UCMH lensing event would provide 
a measurement of the direction and speed of the UCMH's 
motion on the sky. As discussed in Ref. (88|, this infor- 
mation could be used to predict the lensing of stars along 
the UCMH's projected trajectory. Subsequent observa- 
tions of these stars using ground-based telescopes could 
then provide confirmation that an UCMH was responsi- 
ble for the detected event. 



D. Lensing Probabilities 

To determine the likelihood of observing a lensing 
event, one must first calculate the total lensed solid angle 
in the sky A to tai- This is equal to the sum of the indi- 
vidual Aionsod over all UCMHs, if they are rare and small 
enough to assume that they do not overlap. If UCMHs 
originate from a strongly localized enhancement in the 
primordial power spectrum, they have the same initial 
mass Mi. Taking their growth factors g to be equal, 
we can extend this assumption to the current mass Mh- 
Assuming UCMHs are homogeneously distributed with 
respect to the local dark matter density (pdm — 0.4 GeV 
cm -3 ), their local number density is n U cmh = pdm/Mh- 
If the fraction of dark matter in UCMHs today (/o) is 1 
(occurring for / cq ^ 1/300), then n ucmh = f eq Pdm/Mi. 
Otherwise, n ucmh = f Pdm/Mh ^ /o/?dm/(300M;). We 
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FIG. 3: Lensing probabilities as a function of initial UCMH mass Mi, assuming that / cq = 1 and the observed stars have an 
average distance d s — 2kpc. As in FigfJl the four lines correspond to different values of S m i n , and the three panels show the 
constraints for three different UCMH density profiles. 



then have 



Aotal = 47T n ucmh / Sensed (di ) d[ d(di) . (16) 
Jo 

The probability that any given star in the sky is be- 
ing lensed by a UCMH is P cnsc d = ^totai/Asky, where 
A s ky — 5 x 10 11 as 2 is the total area of the sky. Note 
that A\ cnsc d{di) depends on di only through the factor 
(1 — di/d s )d^ 2 , so Piensed is linear in d B . In light of 
this relationship, we adopt an averaged source distance 
of d s — 2 kpc, which conservatively represents the full 
Gaia discovery space including the large number of dis- 
tant giant stars at the brightness levels we target. 

Lensing probabilities are shown in Fig. [3] for f cq = 1. 
As in Fig. [2J the probability drops with larger cores and 
larger S mm , and there is a minimum required halo mass 
Mi for there to be any probability of a star being lensed. 
The probability rises sharply as M t increases beyond this 
minimum, again reflecting the behavior of the individual 
cross-sections. However, because of the n ucm h factor, the 
probability falls for larger halos as their increasing rar- 
ity begins to dominate over their larger individual cross- 
sections. The result is that there is an initial mass Mj 
around which a given detection scenario is most sensitive. 



IV. CONSTRAINTS ON THE PRIMORDIAL 
CURVATURE PERTURBATION 

The expectation value for the number of observed lens- 
ing events is A^bs x Piensed, where N QDS is the number of 
stars sampled. If a single UCMH lensing event is detected 
in such a sample, we can invert the Poisson cumulative 
distribution function to obtain a 95% confidence lower 
bound on the lensing probability: Power = 0.355/iVobs- 
Conversely, the absence of UCMH lensing events implies 
a 95% confidence upper bound on the lensing probability: 
Pipper = 2.996/A' r bs. To determine a feasible value for 
N bs, we start with the Gaia astrometric performance as 



estimated at the time of the Gaia Mission Critical Design 
Review. 3 

The Gaia astrometric performance is approximately 
constant for stars brighter than G=12 (where G is the 
stellar magnitude in the Gaia passband) , and so there is 
no penalty for including all stars down to that bright- 
ness. This cut-off leads to a target list of at least several 
million stars, based on the G-band all-sky star counts 
detailed in Gaia Technical Note GAIA_ML_022 and Ref. 
|l!2j | . To ensure a low number of false positives with that 
number of targets, we require a 6a detection. In this sce- 
nario, our canonical value of S'min — 256 fi&s corresponds 
to a single-measurement precision of 29 /ias. This perfor- 
mance is reached by Gaia for a stellar brightness of ap- 
proximately G=12.5, for which the G-band all-sky star 
counts predict approximately 7 x 10 6 target stars. We 
expect some targets to be rejected during the calibration 
phase because of companions and other effects, so we re- 
duce the estimated AT bs to 5 x 10 6 . With this sample 
size, the 6cr limit for S'min implies 98% confidence in a 
single detection. 

In order to translate constraints on the lensing prob- 
ability into constraints on the initial UCMH mass frac- 
tion, we note that P on sed has a twofold dependence on 
/cq if fcq 1/300. The first is a simple linear depen- 
dence originating from the number density n ucm h oc / eq . 
The second is a more complex term originating from the 
growth of the truncation radius r t as the initial UCMH 
mass fraction decreases, which increases the post-equality 
growth factor g = l// eq . This leads to a weaker restric- 
tion on the lensing cross section j4i onS cd, which must fall 
completely inside the truncation radius. However, for the 
ranges of S m m and Mi under consideration, the lensing 
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FIG. 4: The upper bound on the fraction of dark matter in UCMHs at matter-radiation equality resulting from the nondetection 
of UCMH lensing events in a survey of 5 million stars. The solid lines show the upper bound on / eq given a minimum detectable 
signal of 256 /ias, which corresponds to a 6<r detection by Gaia. The dashed and dotted lines show the upper bounds for more 
sensitive surveys, with Smin = 64 /ias and S'min = 16 //as, respectively. As in Fig[2] the three panels show the constraints for 
three different UCMH density profiles. 



signal falls off so sharply with increasing star-lens sepa- 
ration that this latter factor is negligible. We can thus 

Safely aSSUme that / eq ~ Plcnscd(/eq)/-Plensed(/cq = !)• 

If iVobs stars are monitored and no UCMHs are ob- 
served, then we can place an upper bound on / cq : / cq < 
(2.996/iV bs)/JWed(/eq = !)• These upper bounds on 
/ C q are shown in Fig. [4] We see that a survey of 5 
million stars with S mnl — 256/ias can constrain / eq 
to be less than 0.5 for UCMHs with masses between 
1M Q <, M t <, 20M Q if (av) < 3 x IO^ch^s" 1 . Fig- 
ure [4] also shows how future surveys with higher astro- 
metric precisions could improve this constraint. Even 
with S min = 16 /ias, however, we can only constrain 

/ eq £ 0.007. 

The / eq constraints shown in Fig. @] assume that the 
inner density profiles of UCMHs are not significantly dis- 
turbed from the time of their collapse to the present 
day. Since UCMHs have high central densities with steep 
density profiles and small constant-density cores, it is 
highl y likely th at they survive accretion by larger halos 
HBlmGil- However, if / eq £ 0.01, then UCMHs do 
not grow in isolation before falling into larger, more dif- 
fuse halos. Instead, UCMHs would interact with other 
UCMHs, and the outcome of such interactions has not 
been investigated. A compete analysis of the survival 
probability for UCMHs with / eq <; 0.01 lies beyond the 
scope of this work, but we note two reasons to expect 
that such an analysis will not affect our primary con- 
clusions. First, UCMHs with Mj ~ 5M Q only generate 
astrometric microlensing signals greater than 256 /ias if 
Pydi & 5 x 10~ 4 pc. Given the UCHM's steep density 
profile, most of the mass beyond this radius could be 
stripped from the UCMH without significantly changing 
its astrometric microlensing signature. Second, we can 
calculate how UCMH destruction affects our bounds on 
the primordial power spectrum. If s is the fraction of 
UCMHs that survive to the present day, then the bounds 
on / eq shown in Fig. 3] are bounds on the product s/ eq . 



We will show below that a Gaia-like survey is capable of 
significantly constraining the amplitude of the primordial 
power spectrum as long as s <; 0.3. 

We assume that the initial density perturbations are 
Gaussian and use the Press-Schechter formalism jl20l | 
to translate the constraints on the initial UCMH mass 
fraction / eq into constraints on the mean-squared ampli- 
tude of dark-matter density fluctuations within a sphere 
of radius R at horizon crossing: a^ OI (R). If we assume 
that the UCMHs do not accrete prior to matter-radiation 
equality, then / oq equals the fraction of the dark mat- 
ter contained in regions with dark-matter overdensities 
greater than the minimum matter overdensity required to 
form an UCMH (<J m in) and smaller than the matter over- 
density required to form a primordial black hole (~ 1/4, 
corresponding to a radiation overdensity of 1/3) [59(: 



/e q (M/ 



,1/4 


r s 2 - 




L, Exp 




dS, 



( 17 ) 

where R is the comoving radius containing a dark mat- 
ter mass Mi [see Eq. ([2])]. We are primarily interested in 
values of a^ OI that are much less than 1/4, so / cq is insen- 
sitive to the upper limit on this integration. In contrast, 
/ q depends very strongly on the value of 5 mm . 

Several studies of UCMHs assume that 5 m ; n ~ 10~ 3 
for all UCMH masses [1, H, @, but this value for S m i n 
is only correct to an order of magnitude, and it ignores 
the scale-dependence of <5 m i n . Since smaller-scale per- 
turbations enter the Hubble horizon prior to larger-scale 
perturbations, they have more time to grow and reach 
the collapse threshold before a given collapse redshift. 
Therefore, if we define S m i n (R) to be the overdensity at 
horizon crossing required for the region to collapse prior 
to a given redshift, <5 min will decrease slightly as R de- 
creases. A precise calculation of <5 m i n as a function of scale 
and collapse redshift was recently presented in Ref. H, 
and we briefly summarize their result in the Appendix. 
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FIG. 5: The bounds on the primordial perturbation amplitude resulting from a survey of 5 million stars. The left axis gives 
the amplitude of the primordial curvature power spectrum, while the right axis gives the root-mean-squared amplitude of dark 
matter density fluctuations within a comoving sphere of radius R = fc -1 at horizon entry. The solid curves show the upper 
bounds that would follow from the failure to detect lensing events at different values of Smi n . The dotted curves shows the 
lower bound that would result if an UCMH lensing event were detected by a survey with Smin = 256 /ias. The three panels are 
the same as the three panels in Fig. [2] 



Given <5 m i n (i?) and the upper bounds on / eq (Afj) shown 
in Fig. |4l we can use Eq. (1171) to derive upper bounds on 
a^ m (R). These upper bounds are shown in Fig. [5] (see 
right axis). We also show the lower bound on a^ m that 
would result from the detection of one UCMH lensing 
event by Gaia with S m in = 256 /ias. 

The translation of constraints on a^ OI (R) to con- 
straints on the primordial power spectrum Vn(k) is 
model-dependent; it depends on the scale-dependence 
of Viz{k)- We assume that Vjz is locally scale- 
invariant (i.e. that it does not vary significantly when 
k changes by a factor of a few), and we constrain its 
amplitude at different scales. With this assumption, 
al m (R) ~ 0.908P-R (k = R' 1 ) @; a brief derivation of 
this relation is given in the Appendix. 

The left axis of Fig. [5] shows the upper bounds on 
V-jz^k) that follow from the absence of UCMH lensing 
signals in a survey of 5 million stars. A Gaia-like survey 
with 5 m i n = 256 /ias probes the primordial power spec- 
trum in a fairly narrow band around k ~ 3000 Mpc . If 
the core radius of the UCMH is sufficiently small, Gaia 
could set an upper bound Vtz & 1Q~ 5 at these scales. 
Fig. [5] also shows how surveys with higher astrometric 
precision could further strengthen this constraint over a 
wider range of scales. Alternatively, if Gaia observes a 
lensing signal from an UCMH, we can put a lower bound 
on Vtz] this lower bound also is shown in Fig. [5j 

As illustrated in Fig. [5j the Vtz constraints derived 
from astrometric microlensing depend on the properties 
of the dark matter particle; the upper bounds on Viz in- 
crease as (av) /m x increases because the lensing signal is 
smaller for UCMHs with larger core radii. The only con- 
straint on Vtz for k > 100 Mpc -1 that does not depend 
on the properties of the dark matter particle is the bound 
from the null detection of PBHs: Vn & 0.05 at k ~ 3000 
Mpc~* [6l|. For small UCMH cores, the constraint from 
Gaia (Smin = 256 /ias) is over three orders of magnitude 
stronger! 



We now revisit the issue of UCMH survival. As long 
as s/eq 0.1, decreasing s has little impact on the up- 
per bounds on Viz(k) shown in Fig. [5J which were cal- 
culated assuming that all UCMHs survive (s — 1). For 
example, if s = 1, no UCMH detections in a survey of 
5 million with S m i a = 16 /ias implies Vtz < 9.5 x 10~ 7 
for fc = 1.2 x 10 4 Mpc -1 ; if s = 0.1, this upper bound 
increases to Vtz < 2.1 x 10~ 6 . The constraints are more 
sensitive to changes in s for larger values of s/ eq ; setting 
s = 0.1 for a survey with S m - m = 64 /^as increases the up- 
per bound on Vjz at fc = 6200 Mpc^ 1 by an order of mag- 
nitude. For a Gaia-like survey with ^min = 256 /ias, set- 
ting s = 0.3 implies Vn < 6.3 x 10~ 4 for fc = 3500 Mpc" 1 . 
While weaker than the upper bound shown in Fig. [5j this 
upper bound on Vtz is still significantly stronger than the 
constraint derived from the absence of PBHs. 



V. UCMHS AS GAMMA-RAY SOURCES 

While they are a significant improvement over the 
constraints from PBHs, the upper bounds on VTz(k) 
shown in Fig. [5] are all weaker than the bounds de- 
rived from the fact that Fermi-LAT satellite has not 
observed gamma-rays from dark matter annihilation 
inside UCMHs, but these constraints depend on the 
properties of the dark matter particle. Ref. Q as- 
sumed that m x = 5 TeV, (av) = 3 x 10 -26 cm 3 s _1 , 
and that all the dark matter particles annihilate to 
bb pairs. With these assumptions, they showed that 
the most stringent bounds on the UCMH population in 
the mass range 2 x 1O _6 M <; Mj <, 3 x 1O 5 M follow 
from Fermi-LAT's failure to detect individual UCMHs 
within our galaxy as gamma-ray point sources. The re- 
sulting constraint on the primordial power spectrum is 
Vtz ^ 3 x 10~ 7 for fc ~ 3000 Mpc" 1 . Constraints on 
Vtz from astrometric microlensing by UCMHs probe a 
different region of dark matter parameter-space; as the 
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FIG. 6: The value of (av) below which constraints on / oq from astrometric microlensing are more stringent than the constraints 
from gamma-ray emission. The three panels show different values of S m in, and in all cases, we assume that the astrometric 
microlensing constraints are derived from a survey of 5 million stars that did not detect lensing events from UCMHs. The three 
curves in each panel show the results for different dark matter particle masses. 



ratio (av)/m x decreases, the Fermi-LAT constraints get 
weaker, while the constraints from astrometric microlens- 
ing either get stronger or are unaffected. In this sec- 
tion, we find the region of dark matter parameter-space 
in which the astrometric microlensing constraints on Viz 
surpass the constraints from Fermi-LAT. We will see that 
the constraints shown in the left panel of Fig. [5] can be 
considered conservative; if the dark matter annihilation 
leads to larger UCMH core radii, then Fermi-LAT places 
a stronger upper bound on Vu- 

An UCMH located at a distance d produces an ob- 
served gamma-ray flux 



T = 



l * dN 
, dE' 



dE 



rV0)dr , (18) 



where i?th is the minimum energy required for detec- 
tion, r t is the radius of the UCMH [see Eq. and we 
have assumed that the annihilation proceeds via a single 
channel with a differential photon yield given by dN /dE. 
For our assumed UCMH density profile [Eq. (|12p] with 
n > r c , 



rt 1fi 

2 2 / \ i u 2 S 

r p (r)dr ~ j^:Po r c- 



(19) 



Since r t 3> r c for all UCMHs, we see that the flux only 
depends the UCMH's core radius and density. If the 
core radius of the UCMH is set by dark matter annihi- 
lation, then po depends only on the properties of the 

1/3 

dark matter particle, but r c oc M i ; it follows that 
pfrc oc Mi(m x / (<jv)) 2 / 3 . Conversely, if the annihilation 
core is smaller than the core produced by nonradial infall, 
then p\r\ oc M lm . 

Following Ref. Q, we define d oos to be the distance at 
which the gamma-ray flux from an UCMH equals 4 x 10 -9 
photons cnr 2 s _1 with _E t h = 100 MeV; this is the flux 
required for a 5<r detection of a point-source after one 
year of observations by Fermi-LAT. If all of the dark 
matter is presently contained in UCMHs (/ eq 1/300), 



then the expected number of UCMHs that Fermi-LAT 
should detect is f cq M d<dobs /M u where M d<dobs is the 
mass of the dark matter enclosed in a sphere with ra- 
dius dobs centered at our location. Since Fermi-LAT 
has not detected any UCMHs after one year of oper- 
ation, we can conclude that f eq < 2.996Mi/ M d<dohB at 
95% confidence. If the constraint derived from astro- 
metric microlensing is / eq < icq. am, then the astromet- 
ric microlensing constraints are more powerful than the 
constraints from Fermi-LAT if 



C?obs < 



/ 3 x 2.996-Mi 



1/3 



(20) 



V 47T/eq,AMPdm , 

All the upper bounds on / eq shown in Fig. 0] require 
dobs < 33 pc to surpass the Fermi constraints, so our as- 
sumption that the dark matter density is uniform within 
a sphere with radius d b s is justified. 

Since <i bs depends on the properties of the dark mat- 
ter particle, the upper bound on d bs given by Eq. (|20)) 
implies an upper limit on (av) for a given value of m x . 
These upper bounds on (av) are shown in Fig. [5] To de- 
rive these limits, we used values for d n hs( Mj) ca lculated 
using an extended version of DarkSUSY |5l.ll2lj. assum- 
ing all the dark matter particles annihilate to bb pairs, for 
three values of (av) and m x 122 1. 4 Since d Q bs J 71 / 2 , 



dobs oc (av) 1 / 6 if the UCMH core is set by dark matter 
annihilation, and d bs oc (av) 1 / 2 if the UCMH core is set 
by the nonradial infall; we used these scalings to obtain 
d hs(Mi, (av)) functions for each value of m x . We then 
found the value of (av) that saturates Eq. (|20[) given the 
upper bounds / cq ,AM(-^fi) shown in Fig.[4j 



4 These values for d b s were calculated using a different UCMH 
profile: p(r) = po for r < r c and Eq. for r > r c . With this 
density profile, the integral given by Eq. H19II equals Po r c- To 
obtain d b s for our UCMH density profile [Eq. JT2J], we reduce 
these d b s values by a factor of ^/16/105 before using them in 
our analysis. 
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Figure [6] shows that the dark matter particle must be 
weakly self-annihilating, with (av) <3x 10~ 26 cm 3 s -1 , 
if searches for astrometric microlensing by UCMHs are 
to yield stronger constraints on / eq than Fermi-LAT. For 
these small values of (av) /m x , the core radius set by dark 
matter annihilation is smaller than the core radius set by 
the breakdown of radial infall. Therefore, as long as (av) 
and m x remain unknown, we can safely assume that the 
UCMH density profile is not affected by dark matter an- 
nihilation when calculating the bounds on / eq and V-ji 
that would result from the failure of an astrometric sur- 
vey to observe UCMH lensing events. Such bounds are, 
in effect, the worst-case scenario; although they would 
be weakened if dark matter within UCMH annihilates 
sufficiently to increase the core radius, that rate of dark 
matter annihilation would be sufficient to make UCMHs 
detectable by Fermi-LAT. 



VI. SUMMARY AND DISCUSSION 

UCMHs are small dark matter minihalos that form 
shortly after matter-radiation equality (z <; 1000). Since 
they form earlier than standard minihalos, UCMHs have 
steep density profiles (p oc r -9 / 4 ) and high central densi- 
ties. These properties make UCMHs ideal targets for de- 
tection by astrometric microlensing [88| : when an UCMH 
passes in front of a star, the light from that star is gravita- 
tionally deflected, and the star's image follows a distinc- 
tive trajectory. We calculated the astrometric microlens- 
ing signatures generated by UCMHs, and we found that 
the image trajectories depend strongly on the density 
profile near the center of the UCMHs. We assume that 
UCMHs have a constant-density core; the minimal core 
radius is determined by the breakdown of the radial infall 
approximation that yielded the p oc r~ 9 / 4 density profile, 
but the core radius may be increased if dark-matter an- 
nihilations limit the UCMH central density. Increasing 
the core radius of UCMH decreases its astrometric mi- 
crolensing signature, which implies that astrometric mi- 
crolensing searches for UCMHs are most effective if the 
dark-matter annihilation rate is suppressed. 

The image trajectories produced by microlensing with 
UCMHs are arc-shaped; the star's image moves slowly 
as the center of the UCMH approaches, then it rapidly 
traverses the arc as the center of the UCMH passes by, 
and finally it slows down at the arc's opposite end as the 
UCMH center moves away. Since UCMHs are extended 
objects, it takes several decades for the star's image to re- 
turn to its true position, which makes UCMH microlens- 
ing events easily distinguishable from lensing by point 
masses. The minimal motion of the image prior to the 
UCMH center's passage prompts us to adopt the same 
observing strategy as Ref. [88|: the first two years of 
observations are used to measure the stars' proper mo- 
tions and parallaxes, and then astrometric microlensing 
events over the next four years are detected as deviations 
from the image trajectory predicted from these measure- 



ments. We find that a Gaia-like survey, with a per-epoch 
astrometric precision of 29 microarcseconds for 5 mil- 
lion target stars, is most sensitive to UCMHs with initial 
masses of ~5Mq. Surveys with higher astrometric pre- 
cisions can constrain the abundance of slightly smaller 
UCMHs, but one-microarcsecond precision is required to 
detect astrometric microlensing by UCMHs with masses 
less than ~0.01 M Q . No current or planned astrometric 
instrument can deliver such precision for large numbers 
of targets. Therefore, at least for the near future, astro- 
metric microlensing cannot be used to search for UCMHs 
that result from the enhancement of density perturba- 
tions during QCD phase transition 0, [57j or during a 
matter-dominated era prior to big bang nucleosynthesis 
|58]. Astrometric microlensing can detect UCMHs that 
form from large density fluctuations with wavelengths be- 
tween 0.6 kpc and 6 kpc. 

Since UCMHs form from large-amplitude density fluc- 
tuations (dp/ p *Z 10" 3 at horizon entry) , their abundance 
probes the amplitude of the primordial power spectrum 
on these small scales. If dark matter self- annihilates, the 
high density of dark matter within UCMHs makes them 
bright gamma-ray sources, and they may be detected by 
Fermi-LAT @, @, Q. The fact that Fermi-LAT has not 
detected gamma-ray emission from UCMHs can be used 
to constrain the primordial power spectrum, but this con- 
straint depends on the annihilation rate of dark matter 
particles within UCMHs J8| . The limits on the abundance 
of UCMHs derived from searches for their astrometric mi- 
crolensing signatures would complement the constraints 
derived from the nondetcction of gamma-ray emission 
from UCMHs because the UCMHs produce larger image 
deflections if the dark-matter annihilation rate is low. 
If dark matter annihilation within UCMHs is efficient 
enough to change the UCMH density profile and weaken 
the constraints from astrometric microlensing, then the 
strongest bounds on the power spectrum follow from the 
fact that Fermi-LAT has not yet observed emissions from 
UCMHs . However, if dark matter annihilation does not 
increase the size of the UCMHs' cores, we find that a 
search for UCMHs by Gaia could constrain the ampli- 
tude of the primordial power spectrum to be less than 
10~ 5 on scales of k ~ 3500 Mpc -1 , an improvement of 
three orders of magnitude over the bound derived from 
the absence of PBHs [6l|. Unlike other constraints de- 
rived from UCMHs, this upper bound on the primordial 
power spectrum does not depend on the properties of the 
dark matter particle. 

When calculating the bounds on the primordial power 
spectrum that result from limits on the local UCMH 
abundance, we assumed that the UCMHs' innermost re- 
gions survive to the present day. UCMHs are proba- 
bly not disrupted during their accretion by larger halos 
because t hey have s mall cores with high central densi- 
ties 0, H, Ill8l Ill9| . The outcome of interactions be- 
tween UCMHs is less certain, and the bound on the local 
UCMH abundance derived from astrometric microlens- 
ing searches is not low enough to ensure that such in- 
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teractions are uncommon. If UCMH-UCMH interactions 
increase the UCMHs' core radii or decrease their central 
densities, then our projected bounds on the primordial 
power spectrum are weakened. We encourage further 
study of UCMH survival in scenarios in which UCMH in- 
teractions are common to determine exactly how UCMHs 
are effected by such interactions and what impact this has 
on their astrometric microlensing signatures. We note, 
however, that Gaia could significantly improve the PBH 
constraint on the primordial power spectrum even if up 
to 70% of UCMHs that formed initially cannot generate 
an observable astrometric microlensing signal today. 

Our analysis could also be refined by considering more 
sophisticated techniques for detecting astrometric lensing 
events. Our detection strategy is designed to evaluate 
the level of astrometric signal produced by UCMHs. An 
optimized detection strategy, perhaps based on matched 
filters, would provide greater sensitivity and improved 
false positive rejection. Finally, we note that gravita- 
tional lensing by dark matter minihalos will also generate 
time delays, w hich could be detected using high-precision 
pulsar timing [l23l Il24j ; it would be interesting to inves- 
tigate how pulsar timing may constrain the abundance 
of UCMHs. 

The upper limit on the amplitude of the primordial 
power spectrum on small scales derived from the absence 
of PBHs is a valua ble tool i n our quest to understand in- 
flationary physics p5l4l3l1 . Limits on the UCMH abun- 
dance can provide a much stronger bound on the ampli- 
tude of the primordial power spectrum on small scales. 
Currently, the only constraints on the UCMH popula- 
tion are derived from their possible emission of dark mat- 
ter annihilation products; these constraints do not apply 
if the dark matter annihilation rate is sufficiently sup- 
pressed. We have shown that astrometric microlensing 
offers a promising alternate method of detecting UCMHs 
and thus constraining inflationary physics. 
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Appendix A: From / eq (M 4 ) to Pn{k) 

In this Appendix, we briefly review how bounds on / eq 
translate into bounds on the amplitude of the primordial 
power spectrum P-R.(k). From Eq. (fTT|) . we see that we 
must first calculate the minimum matter overdensity re- 
quired to form an UCMH (<5 m i n ). Once we have <5 m in(-R), 
we can find the value of er^ or (i?) that yields the desired 
value of / eq . If we then assume a specific scale depen- 
dence for Pjz(k), we can use <j 2 Q1 .(R) to obtain P-n(k). 
This procedure is described in detail in Ref. [8j for several 
models of Pn(k). We briefly summarize the key results 
here. 

During matter domination, a region collapses when 
the linear overdensity reaches the collapse overdensity: 
Scoll = (3/5)(3tt/2) 2 / 3 ~ 1.686. The amplitude of 
subhorizon density perturbations during the matter- 
dominated era is proportional to the amplitude of the 
curvature perturbation 1Z at horizon entry: 



5 x (k,z) 



2 k 2 



5H 2 tt 



M 



-T(k)V(z)U(k) 



(Al) 



where T(k) is the transfer function normalized to unity 
on large scales, and 



V{z) = 



fl + z Rn ) 



1 



l + z 



- 2 



l + z 



is the grow th function after matter-radiation equality 
(e.g. jl32| ). We can use Eq. (|Al|) to find the curva- 
ture perturbation required to make 5 x (k, z co \\)_= <5 co n at 
some collapse redshift z C oii- We follow Ref. Q and set 
«boU = 1000. 

Next, we need to calculate the dark matter density 
perturbation at horizon entry that corresponds to this 
curvature perturbation; this is <5 m i n . The overdensity at 
horizon entry is a gauge-dependent quantity, however. 
Since 5 co \i was evaluated in the density rest frame (total 
matter gauge), and the relevant modes for UCMHs en- 
tered the horizon during radiation domination, Ref. Q 
evaluates <5 m i n in the rest frame of the radiation. In this 
gauge, the dark matter density perturbation during ra- 
diation domination is given by 



5 x (k,a) = 6K{k) 



lnfl + 7£ -i-Ci(tf) + ^ 



(A3) 



where is the Euler-Mascheroni constant, Ci is the co- 
sine integral function, and 9 — k/{\fiaH). Evaluating 
this expression at horizon entry (k = aH) and setting 1Z 
to the minimum value required to make S x (k, 2 C oll) = ^coll 
yields 



[ 'min (k, Zcoll) 



5 H 2 fl M 0.988 

6 c ° n k 2 T(k)D(z co n)' 



(A4) 



For k = 5.1 x 10 4 Mpc 1 , which corresponds to Mj = 
10 _3 Af Q , (5 min = 2.2 x 10~ 3 . As expected, 6 m i n is slightly 
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larger for larger UCMHs; for M t = 1OOOM , S min = 
3.4 x 1CT 3 . 

The primordial power spectrum determines a 2 {R,t): 

r°° d/fc 
<j 2 (R,t) = J F 2 (kR)V x (k,t) T , 

where F(x) = 3x~ 3 (since — cccosx) is the Fourier trans- 
form of the top-hat window function, and T 2 (9) = 5 x /9 2 
with 5 X given by Eq. (|A3|) . If we evaluate a 2 (R,t) at 
horizon crossing [R = (aH)^ 1 ] and change the integra- 



tion variable to x = kR, we have 



<AR) = l jf° ^T 2 X (e = F 2 {x)V n (k = |) dx 

(A6) 

The factor x 3 T 2 F 2 serves as a modified window function: 
it is peaked at kR ~ 2.4, and if V-R.(k) is scale-invariant, 
80% of the integral's value comes from integrating the 
range 1 < kR < 4. We assume that V-r is nearly constant 
within this limited range of scales. In this case, V-jz(k) 
can be pulled outside the integral, and we are left with 
al OT (R)~Q.9Q8V 1 i(k = R- 1 ). 
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